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Abstract. 

Phenomena such as pulsar frequency glitches are believed to be attributable to 
differential rotation of a current of "free" superfluid neutrons at densities above the 
"drip" threshold in the ionic crust of a neutron star. Such relative flow is shown to be 
locally describable by adaption of a canonical two fluid treatment that emphasizes the 
role of the momentum covectors constructed by differentiation of action with respect 
to the currents, with allowance for stratification whereby the ionic number current may 
be conserved even when the ionic charge number Z is altered by beta processes. It 
is demonstrated that the gauge freedom to make different choices of the chemical basis 
determining which neutrons are counted as "free" does not affect their "superfluid" 
momentum covector, which must locally have the form of a gradient (though it does 
affect the "normal" momentum covector characterising the protons and those neutrons 
that are considered to be "confined" in the nuclei). It is shown how the effect of 
"entrainment" (whereby the momentum directions deviate from those of the currents) 
is controlled by the (gauge independent) mobility coefficient JC , estimated in recent 
microscopical quantum mechanical investigations, which suggest that the corresponding 
(gauge dependent) "effective mass" 777.* of the free neutrons can become very large in 
some layers. The relation between this treatment of the crust layers and related work 
(using different definitions of "effective mass") intended for the deeper core layers is 
discussed. 



1 Introduction 



As a prerequisite for the quantitative analysis of the role of differential rotation in the 
angular momentum transfer mechanism that is thought [I] to be responsible for phe- 
nomena such as pulsar glitches, the purpose of this article is to describe the adaptation 
to the special circumstances pertaining to the inner crust of a neutron star of the kind 
of non relativistic two constituent fluid formalism that has already been applied [21 E] 
to the global description of the main bulk of the star, in which the basic constituents 
are relatively moving currents of protons and neutrons, using a treatment of the non- 
dissipative kind [3J [5] that has recently been developed in a variational formulation 
jHllZllH] that, although non-relativistic, can nevertheless be very conveniently expressed 
in a fully covariant manner (A non-dissipative treatment is useful as a good first 

approximation in view of the superfluidity of the neutrons, but in any case, after the 
appropriate action function has been obtained in the manner described below, it will 
be possible to use it as an equation of state in a less idealised model of the kind [TTj 
needed to allow for dissipation and for the relativistic effects JE] that are important at 
a global level.) 

Such a two fluid treatment does not take account of the third independent current, 
namely that of the electrons, that is of dominant importance for electromagnetic effects, 
but which makes a relatively negligible contribution to the Newtonian mass transport 
effects under consideration here. Allowance will however be made here for an important 
special feature of the crust layers (providing stratification that contributes to stability 
against convection) namely the clustering of the protons into ions, whose number current 
can remain conserved even after allowance for possible variation of the ionic charge 
number Z> due to weak interactions whereby protons are transformed into neutron or 
vice versa. Another special stabilising feature of the crust that will also need to be 
allowed for, but that will be left for subsequent work is that of the (relatively small) 
stress anisotropy that can arise from the elastic solidity jTJj of the crust, or from strong 
magnetic fields [T3| . 

The present article is particularly concerned rather with allowance for the effect of 
relative "entrainment" fH] between on one hand the effectively free neutrons that will 
be present above the "neutron drip" density threshold (of the order of 10 11 g/cm 3 ) and 
that are expected to behave as a superfluid except at very high temperature, and on 
the other hand the protons, which (together with the ions to which they are confined) 
will behave as an effectively "normal" fluid in the crust, though they may behave as 
another superfluid in the deeper layers at densities above the nuclear threshold of the 
order of 10 14 g/cm 3 at which the ions merge. 

One of the main goals of this article is to clarify the relationships between different 
kinds jUE] of definition of "effective mass" that have been introduced by various authors 
as a quantitative measure of the "entrainment" effect, and particularly to provide a 
discussion of the way such masses depend on the gauge freedom involved in the choice 
of a chemical basis for the purpose of specifying just which neutrons are considered to 
be free. In the homogeneous (unclustered) phase above the nuclear density threshold 
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prevailing in the liquid core the neutrons will all be effectively free, while in the outer 
crust below the "drip" threshold it is clear that none of them will be free, in the 
"operational" sense of being able to penetrate the potential barriers between nuclei in 
an astrophysical timescale shorter than the age of the universe. In the inner crust at 
densities just above this threshold, neutrons in low and intermediate energy states will 
still be effectively confined, but there will be a clear cut critical energy above which 
they will be effectively free to travel over ionic separation distances on a microscopically 
short timescale. 

The problem is that, at deeper levels in the crust, such an "operational" discrim- 
ination between "confined" and "free" neutron states will no longer be so clear cut, 
because the ionic wells will get too close for the exponential suppression outside to be 
fully effective, so that there will be marginally bound states with intermediate penetra- 
tion timescales that are macroscopically long but cosmologicaly short. This means that 
there will be some degree of ambiguity that needs to be resolved by a more or less arbi- 
trary choice of a chemical gauge whereby the baryon number A of the ionic clusters 
is prescribed as a function of their (depth dependent) charge number Z . Fortunately, 
as will be shown below, the most important quantities in the two fluid formalism at 
a local (intervortex) level notably the action density A with its concomitant energy 
function, and in particular the (irrotational) superfluid neutron momentum, all turn out 
to be robustly independent of which particular ( "operational" or other) gauge criterion 
is used. 

Following this discussion of the different ways of formulating the generic treatment 
of entrainment in the two fluid formalism, a more particular purpose of this work is 
to show how quantitative values of the relevant "effective mass" are provided by the 
corresponding (gauge independent) "mobility coefficient" /C for which quantitative 
estimates are obtainable from the kind of microscopic analysis described in our im- 
mediately preceding work [16J. This analysis is based on a non-relativistic quantum 
mechanical formalism of the type introduced by Oyamatsu and Yamada in a local 
frame with respect to which the ionic crust lattice is at rest, so that it gives rise to a 
static effective potential having a periodic dependence on the cartesian space coordi- 
nates x 1 . The analysis is carried out in terms of single particle wave functions subject 
to Bloch type boundary conditions (meaning that they deviate from ordinary period- 
icity by a phase factor of the form exp{i/CjX*} ) of the kind that is commonly used 
for electrons in ordinary solid state physics, but that has not yet been systematically 
employed for analogous problems involving neutrons. In a follow up article ^Hj we have 
shown how the method can be extended to include the moderate adjustments arising 
from the effect of BCS type superfluid pairing. 

The mobility coefficient /C , that is provided by this analysis, has a value which is 
interpretable as the ratio of the density n of "free" ("conduction") neutrons to their 
effective mass m+ ■ It was suggested by our preliminary results fH] and has been 
confirmed by more recent work that this value is likely to be much larger than had 
been expected. Even when measured with respect to an "operational" gauge as shown 
in Figure 0] in which the number density 71 of neutrons that are counted as "free" 
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allows only for those that are not bound inside ionic nuclei, the quantity is likely 
to become extremely large compared with the ordinary neutron mass m in the middle 
layers of the inner crust. It would of course be even larger with respect to the very 
simple "comprehensive" gauge in which all the neutrons are counted as "free" so that, 
as shown in Figure El the corresponding "effective mass" would in any case diverge as 
the density decreases to the "neutron drip" value. 

The purely Newtonian (i.e. non relativistic) framework within which the work of 
the present article is carried out should be adequate, on a local scale, as a fairly good 
approximation in the density regime under consideration. This lower crust regime 
extends from the "neutron drip" threshold to the transition to a purely fluid (at low 
temperatures actually superffuid) mixture of neutrons and protons at the base of the 
crust where the mean density approaches that of ordinary nuclear matter (at a value 
of the order of about 10 14 g.cm~ 3 .) In a global representation of the neutron star (on 
scales of the order of a kilometer or more) there will be significant deviations from flat 
geometry whose quantitative treatment would require the use of General Relativity, but 
the discussion here will be limited to a local neighbourhood that is sufficiently small (on 
scales of the order of a centimeter or less) to be treated as approximately homogeneous 
and geometrically flat. Even on such a purely local scale it would be necessary to use 
a special relativistic (Minkowski space) description if we were concerned with effects 
involving the electrons, but the rest mass of the neutrons, to which the present analysis is 
restricted, is so much higher that a purely Newtonian analysis should be quite adequate 
in the moderate energy range (up to a few tens of MeV) that is involved. This contrasts 
with the state of affairs in the fluid core of the neutron star, where estimation of the 
relevant "entrainment" between (superffuid) neutrons and protons requires a different 
kind of treatment [2D] i n which allowance for relativity is important, but where the 
effect of the entrainment is much more moderate. An implication is that to improve 
on previous work [21] that ignored the effects both of relativity and of entrainment in 
the analysis of angular momentum transfer in the crust, the effect of relativity is of 
secondary importance, and that the first priority is to include the relevant allowance 
for the effective mass enhancement due to entrainment [22]- 

2 Generic two-fluid models with stratification 
2.1 The canonical master function 

In astrophysical contexts for which a non-relativisitic Newtonian mechanical description 
is adequate it will usually be sufficiently accurate to represent the relevant distribution 
of mass just in terms of a conserved current of baryons characterised by a single constant 
mass parameter, m say, that can be taken to be the standard atomic mass unit or 
simply the proton mass which is different only by a fraction of a per cent and so is near 
enough for practical purpose. It is convenient for many purposes to use a 4-dimensional 
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notation scheme in which, for example, the total baryon number density 77-b 

say, and the corresponding baryon current density components with respect to space 
coordinates X 1 (i = 1, 2, 3) are combined to form a 4-current [y = 0, 1, 2, 3) 
whose time component is given by n-^ = 77-b ■ This means that the total baryon 
conservation law, 

dn h /dt + Vini = 0, (2.1) 
will be expressible more concisely as 

V„< = 0. (2.2) 

(Throughout this paper, summation over repeated indices is assumed). 

Our purpose here is to consider situations involving relative motion between an un- 
charged contribution that will be attributable mainly to neutrons (but perhaps 
including other neutral hyperons such as the A in the inner core regions) and an elec- 
trically charged contribution that will similarly be attributable mainly to protons. 
(It is to be understood that overall electrical charge balance is ensured by a non bary- 
onic current attributable mainly to electrons, whose mass contribution is so small that 
it can be neglected.) The two independent baryonic contributions combine to give the 
conserved total baryon current 

< = < + <• (2-3) 

In dynamic processes on short timescales the neutron and proton currents may each 
be considered to be separately conserved, meaning that the 4-divergences Vj/Ti^ and 
Vj/Tlp will both vanish individually, but in long term "secular" evolution processes it 
will be necessary to allow for the possibility of converting neutrons to protons or vice 
versa by weak interaction processes. 

The present discussion will be concerned with cases in which the evolution of the 
separate contributions and can be described by a multiconstituent fluid model 
of the kind governed by a master function A ma t , that acts as the material (meaning non 
gravitational) part of a Lagrangian in effectively conservative applications [HJ E2] , but 
that is also compatible with allowance for potentially relevant dissipative effects [TTj . 
Such a multifluid description does however entail neglect of the extra energy contribu- 
tions due to possible elastic solid deformations [Hj, and to frozen in magnetic fields [TH] . 
whose inclusion would require a more elaborate treatment. As well as depending on 
the separate currents and it is useful to allow for the potentially important 
effect of stratification due to baryonic clustering forming ionic nuclei characterised by 
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a charge number Z and hence by a number density that is expressible as 

m = n p /Z, (2.4) 

wherever the density is below the saturation density of the order of 10 14 g.cm 3 . This 
means that the generic variation of the master function will be given in terms of corre- 
sponding partial derivative coefficients by an expression of the form 

£A ma t = -/J 1 Sni + //1 5n£ + yP v 6n^ , (2.5) 

or equivalently in the less compact form of a traditional 3+1 space time decomposition 
by an expression of the form 

£A mat = -/i 1 Sni - /i n 8n n - ^ p Sn p + 5n l n + ^ Sn^ , (2.6) 

in which flf and /i? will be respectively interpretable as the mean momentum per 
particle of the neutrons and of the protons, while the quantities 

M n = -fi , ^ = , (2-7) 

will be interpretable as the corresponding neutronic and protonic chemical potentials, 
while the coefficient \l is an ionic cluster potential, whose gradient, if any, represents 
the effect of stratification. This potential can be used to construct a corresponding 
4-momentum covector that is given in terms of the gradient of the Newtonian time 
coordinate X° = t by the formula 

fi v = -H l t Vl t v = V v t = 5l, (2.8) 

which means that it has vanishing space components, (I a == 0, as an expression of 
the fact that (in the Newtonian limit) the clustering contributes only to energy but not 
to mass. 

The models we are considering are intended for application to neutron star matter 
after it has cooled down sufficiently for its dynamics to be hardly affected by thermal 
effects, which are not included. These models are to be interpreted as being governed 
by dynamical equations that are formally identical to those of the heat convecting 
thermal model presented in the appendix of the work referred to above, with the 
ionic cluster potential /i 1 taking the place of the temperature , and with the ionic 
number density tl\ taking the place of the entropy density 6 , so that the entropy 
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current 4-vector 6^ is to be replaced here by the ionic number current 4-vector Tlx 
as given in terms of the flow 4-vector of the crust by 

n{ = m < , < = </n p . (2.9) 



In terms of the generalised pressure function 

* = Amat + ni yt - n^l - , (2.10) 

and the Newtonian gravitational potential (J) the corresponding material energy mo- 
mentum 4-tensor ^0] will therefore be given by 

T» v = </xi + n p >P + n>» + W v - mn»t v . (2.11) 

This means that, if the only external force is gravitational, the total force balance 
equation will then be expressible as 

= -mn h V v (f> , (2.12) 

while the condition of conservation of entropy in the conservative case is to be replaced 
here by the ionic cluster conservation condition 

V„< = 0, (2.13) 

which can be expected to be valid, not just in short timescale dynamic processes in 
which there would be separate conservation of the neutrons 

V,< = (2.14) 

(and hence also, according to ([2.2|h of the protons) but even in processes having the 
much longer timescales needed for weak interactions to adjust the relative numbers 
of protons and neutrons so as to achieve the condition of "beta equilibrium", namely 
vanishing of the affinity 

A=(j$-iJi)u?, (2.15) 

as specified with respect to the "normal" rest frame specified by , which is that of 
the protons (and therefore also of the ions) as opposed to that of the neutrons, whose 
superfluidity at low temperatures allows them to retain a relative motion (at densities 
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above the neutron drip threshold) even in a state of exact non dissipative equilibrium. 
On timescales even longer than needed for the beta equilibrium condition, 

.4 = 0, (2.16) 

it is ultimately to be expected that barrier tunnelling processes involving nuclear 
fission and recombination would lead to violation of (j2.13|) in such a way as to efface 
the stratification by allowing the ionic cluster potential fi to tend towards zero, but 
in practice the stratification will commonly survive long enough to have important 
stabilising consequences, particularly in the outer layers of the neutron star. 

To provide a complete system of evolution equations for the dynamic variables 7l\, 
11^, and Tip 1 ', the system consisting of (J2.12j) and (J2.13j) . together with one of the 
alternative possibilities ()2.14|) or (|2.16|) . does not suffice, but must be supplemented by a 
further condition which can simply be taken to be the mesoscopic neutron superfluidity 
condition to the effect that the total neutron momentum covector, 

nl = fil-<f)mU, (2.17) 

should be proportional to a phase gradient, taking the form 

< = W^ n , (2.18) 

in which the phase variable (p n has period TT (not 2tt ) due to the fermionic half 
integer spin of the neutrons JTHJ - (This corresponds to a period having the usual value 
2tt for the corresponding phase if = 2(p n for Cooper pairs with momentum covector 

2< .) 

2.2 Chemical invariance of superfluid momentum 

In order to obtain a formulation that is consistent with traditional usage in the low 
density regime with density less than the order of 10 11 g.cm 3 for which the neutrons 
will also be effectively confined within the nuclei, it is useful to introduce a dependent 
variable, A c , that is to be interpreted as the total number of effectively "confined" 
baryons per nucleus, and that specifies a corresponding "confined" baryon number 
current , which will be given by 

nZ = A c n{. (2.19) 

The "confined" baryon number A c = n c / n\ will evidently lie somewhere in the range 
Z < A c < At, where A\^ = Tl^/fli is the total baryon number per nucleus, which 
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will of course be the same as A c at densities below the "neutron drip" threshold. As 
a corollary there will also be a corresponding "free" neutron current given by 



< = <-<, (2.20) 

which will be non vanishing at densities above the "neutron drip" threshold. In the 
deeper layers of the ionic "crust" region where the clusters are near to stage of merging 
they may become highly deformed to "spaghetti" or even "lasagna" shaped configu- 
rations E3] for which (as remarked above in the introduction) the meaning of the 
absolute numbers Z and A may become rather hazy, but this need not prevent the 
ratio 

a c = A c /Z = n c /n p (2.21) 

of the "confined" baryon number density to the charged baryon number density from 
remaining definable according to some ansatz that might reasonably be required to 
specify a c as a weakly dependent function of Z in such a way that beyond the 
transition to an unclustered homogeneous phase above ordinary nuclear matter density, 
the corresponding "free" neutron current 71^ will simply become the same as the entire 
neutron current n% , while at lower densities it would be determined by the equation 
of state for "cold catalysed" matter whereby all relevant scalars including Z and A c 
are given (in such a way as to minimise the energy density) as functions of TL\, . 

However that may be, subject to the specification of any reasonable (realistic or 

idealised) equation of state for A c as a function just of Z so that the variation of the 

ratio a c {Z} will be given in terms of its derivative a! c = da c / dZ by an expression 

of the form 

5ac = Za' c - **) , (2.22) 

it will be possible, and for some purposes convenient, to make a corresponding chemical 
basis transformation whereby the new (empirically defined) current variables 71% and 
71% are used in place of the original (more physically fundamental) current variables 
71% and n% as the independent variables of the system. In terms of the new chemical 
basis, as specified by the homogeneous linear transformation 

n% = n% + (1 - a c )n% , n% = a c n% , (2.23) 

the generic variation (|2.5jl of the master function and the corresponding energy mo- 
mentum tensor (j2.11j) will be expressible in the exactly analogous forms 

£A ma t = -A 1 6ni + //„ 5n% + 6n% , (2.24) 
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and 

T\ = ntfp v + nful + + " (f> mngU , (2.25) 

with 

* = Amat + ni li l - n»nl - n\^ v , (2.26) 

in which the adjusted ionic cluster potential jd l and the correspondingly adjusted 
analogue 

A = (ji - nlM , (2.27) 

of the affinity (|2.15j) will be related to their untransformed analogues by 

A = A{a c + Za' c ) (2.28) 

and 

fi I = pL I - Z 2 Aa' c , (2.29) 

while the new "free" and "confined" baryon momentum 4 covectors \± v and \J? V will 
be given in terms of their untransformed analogues /i 1 , fl u , /i^,, by the relations 

III = \i% + (1 - OcX/4, - fi%) - ZAa' c t v (2.30) 

and 

fi% = n{ . (2.31) 

This last result is important: it tells us that whereas the momentum covector fi% 
associated with the "confined" part has a rather complicated dependence on the chem- 
ical gauge specified by the choice of the functional dependence of a c on Z, on the 
other hand the momentum covector fj, associated with the "free" neutron current is 
invariant with respect to the change of gauge. The correponding total 

= fJ? v -<f>mtp, (2.32) 

is just the same as its analogue (|2.17J) in the original formulation, so that the super- 
fluidity condition (|2.18|) will be expressible in terms of the the same phase variable ip 
in exactly the same way as before, meaning that it will be given simply by 

n* v = hV v <p. (2.33) 
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The upshot of this is that if the set of dynamical equations is completed by the 
beta equilibrium condition (|2.1fi|) . which according to (|2.28j) will be expressible in the 
language of the new formulation simply by 

A = 0, (2.34) 

then the entire system will be chemically covariant in the sense of having the same 
form regardless of the choice of the gauge function a c {Z} in the chemical basis 
transformation (|2.23j) . 

On the other hand if (as will be appropriate for cases involving relatively short dy- 
namical timescales) the beta equilibrium condition (J2.16)) is replaced by the separate 
neutron current conservation condition (j2.14j) then the complete system will be chem- 
ically covariant only with respect to transformations in which CL C is simply taken to 
be constant, as for example in what we shall refer to as the "comprehensive gauge", 
which is given simply by a c = 1 so that absolutely all of the neutrons are classified 
as "free" , while a less trivial example is that of what may be described as the "paired 
gauge" , which is given by a c = 2 , and which is interpretable as meaning that the only 
neutrons classified as "confined" are those in the tightly bound states that are directly 
paired with corresponding proton states in the nuclei. However, if we wish to use d c to 
provide a more realisticlly "operational" estimate of the total baryon to proton number 
ratio in the nuclei, it will need to be given a non vanishing derivative CL C , and in that 
case the corresponding formulation of the separate conservation condition (I2.14J1 will 
have the patently gauge dependent form 

V„<= -a' c ny v Z . (2.35) 
2.3 Local energy function 

In order to obtain information about the equation of state giving the functional depen- 
dence of the material action density A ma t on the relevant independent variables by 
comparison with the results of a microphysical analysis, it is convenient to work with the 
corresponding energy function. The complete energy momentum tensor (|2.25jl comes 
from a complete action density A = A ma t + A po t in which, as discussed in the pre- 
viously cited work [TOj, the gravitational potential energy contribution is given simply 
by Ap t = —p<fii with p = mn^ so the corresponding energy U = ~T° o will 
have the form of a sum U = U m& t + Up t involving a gravitational contribution 
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U po t = — ^pot that must be subtracted off to feave the purely material part U ma t 
corresponding just to A ma t . It can thus be seen that this material energy density will 
be given (for an arbitrary choice of the chemical gauge function a c {zv } ) by 

U mat = n\ - n\ni - ld c Q n c - * , (2.36) 

which, by ()2.26|) , is evidently equivalent to taking 

U mat = /4 n\ + [l\ n* - A mat . (2.37) 

So long as the relevant velocities are sufficiently small, it will be possible to decom- 
pose the material Lagrangian and energy density in the form 

A ma t = Ai ns + Ad yn , (2.38) 

and 

U ma ± = U[ ns + t^dyn ■> (2.39) 
in terms of a static internal energy contribution 

£Ans = -Ains, (2.40) 

that will depend only on the relevant ionic and (free and confined) baryon number 
densities namely nj, rz-f , n c , together with a dynamic contribution for which the ve- 
locity dependence is homogeneously quadratic so that as for the purely kinetic part [TU] 
the corresponding energy contribution will be the same as the Lagrangian contribution 
from which it is derived, having the form 

A d yn = U dyn = n\ + [i\ n*) . (2.41) 

This means that it will be given as a function of the current components Til and 
by an expression of the form 

Adyn = \m (/C ff n\ n{ + 2/C fc n\ + /C cc n\ n^) , (2.42) 

in which the coefficients , /C cc , and K} c are functions only of 71t, 7T-f, and n c . 
It follows from the Galilean covariance requirements jHJ in the the decomposition ()2.58|) 
to be discussed in the next section that these coefficients must satisfy the conditions 

/C% + /C fc n c = m , JC cc n c + /C fc n f = m , (2.43) 



12 



which means that only one of these coefficients needs to be specified. It follows from 
(I2.42j) that the momenta defined by (|2.fij) will be given as homogeneous linear functions 
of the currents by 

£ = 7ij (/C ff n( + /C fc ni) , = ll3 (/C fc n{ + /C cc n/) . (2.44) 

These relations will be invertible to give the expressions 

n\ = Y J (/C ff ^ + /Cfc /i j) , n l c = 7 « (/Cfc //* + /C cc /4) , (2.45) 
in which the original coefficients will be related to the new ones by 

i^ffi^cc )k"fc2 ' j^ff j^cc i^fc2 ' /Ccc j^fjj^cc j^-fc2 ' (2.46) 

/L A^ — A^ A^ — A^ A^ A^ — A^ 

In terms of these new coefficients, the covariance conditions (|2.47J) are expressible as 

/Cff + /Cfc = — , /C cc + /Cfc = — . (2.47) 
m m 

The expression (|2.42|) is convertible to the dual form 

^dyn = \l 13 (/Cff /4 + 2/Cfc ^ if j + /Ccc Mi Mj) ■ ( 2 - 4 8) 

The requirement that (in order to ensure that the minimum energy configuration is 
the one in which the currents vanish) the dynamical energy density U dy n in expressions 
fl2.42j) or ()2.48j) should be positive definite , implies that the on-diagonal coefficients 
should all be positive, as well as the restrictions 

/Cff /C cc > (/Cfc) 2 , /C ff /C cc > (/C fc ) 2 , (2.49) 

which can be seen from ()2.46j) to entail that the off-diagonal coefficients K} c and /Cfc 
should have opposite signs. 

2.4 Specification of the entrainment coefficient 

In order to relate the canonical approach described above to a description of the tradi- 
tional kind jU in terms of 3- velocities and densities, the dynamical energy density 
has to be rewritten in terms of the relevant current 3-velocity vectors V* and as 
specifed with respect to chosen Galilean rest frame by setting 

n* = n c v l c , n\ = rifV^ , (2.50) 
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which means that v£ will be physically well defined as the "normal" velocity of the 
ions and the rigorously confined proton current, for which we shall have 

ni = ny c , nl = n v v l c . (2.51) 

whereas v£ depends on the chemical gauge function a c {Z} that is chosen to deter- 
mine which of the neutrons are classified as "free" . The chemical gauge dependent but 
Galilean frame independent relative flow velocity vector 

4 = 4 - < , (2-52) 
can then be used for construction of the relative current 3- vector 

n { = mvl, (2.53) 

which as well as being evidently frame independent, also has the noteworthy property 
of being physically well defined in the sense that (like the neutron momentum covector 

£ ______ 

\l i ) it is unaffected by a change of the gauge parameter a c in ()2.23|) . This can 
be seen from the possibility of rewriting it directly in terms of the unambiguously well 
defined total baryonic number density Tli, = -\- n c , and the unambigously defined 
"normal" 3-velocity vector V c , in the form 

= ni - n h v l c , (2.54) 

which is manifestly independent of any prescription for specifying which neutrons are 
deemed to be "confined" and which are deemed to be "free". (Before continuing it is 
to be remarked that subject to substitution of total lepton number current, n 1 ^ say, 
in place of the total baryon number current 77^ above, a precisely analogous relation 
applies in the context of electron conductivity in ordinary solid state physics.) 

We can now proceed to convert the dynamical energy formula to the generic form 
U dyn = \la (p S v} v> + 2// c v\ vi + /. cc ^ V?) , (2.55) 
in which the required mass density matrix components can be read out as 

/ = /C ff ?? f 2 , p tc = JC ic ni n c , p cc = /C cc n c 2 . (2.56) 



This expression (|2.55J) is to be compared with the alternative arrangement whereby 
the decomposition 

P = Pt + Pc , (2.57) 
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of the total mass density p = 777 77b , determines a corresponding decomposition 

tfdyn= ^kin+^ent, (2.58) 
in which the first term is a kinetic energy contribution of the standard form 

^kin = ^(/W + /W) , (2.59) 

while the second term U en t in (|2.58j) is a frame independent "entrainment" contri- 
bution that can depend only on the magnitude Vf c of the relative velocity ()2.52|) and 
must therefore be given in terms of some "entrainment" mass density pi c by 

Uent = 2 Picnic 2 ■ (2.60) 

Substituting ()2.60|) and ()2.59|) in ()2.58|) , and comparing the result with the preceding 
expression (|2.55j) for U dy n , we see that the density coefficients in the latter will be 
given by 

/ = (H + jo fc , p ic = -p fc , p cc = p c + p fc . (2.61) 

It is to be observed that the values of all the density coefficients involved will of 
course depend on the choice of the chemical gauge parameter a c that is used for the 
specification 

Pi = mrii = m(nb — a c n-p) p c = mn c = ma c n p . (2.62) 

of the relevant "bare" mass densities in the decomposition (|2.57jl . 

The foregoing relations can be translated into the terminology of effective masses 
by introducing the values 777^. and 777^ according to the specifications 

p s = m ml , p cc = 77 c , (2.63) 

so that one obtains 

7 Vj = rriivi - vl) + mvl , f'p) = mjfo' - ^) + 7777;/ . (2.64) 

The respective deviations, 777^ c and 777 C £ say, of these effective masses from the 
ordinary baryonic mass 777 will be given in terms of the entrainment density by the 
familiar expressions [Sj 

f* f* Pi. Pi 

777 = 777, — 777 = — - , 777 C f = 777^ — 777 = — - . (2.65) 

77f 77 c 
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It can be seen from (J2.65|) that both effective masses are either larger or smaller 
than m, but one effective mass cannot be larger than m and the other smaller. 
In particular whenever one of the effective masses coincides with the ordinary baryon 
mass, this entails that the other one is also equal to the ordinary mass. 

It is important to distinguish the true (current) 3-velocities used above from the 

pseudo-velocities that are referred to in many older and even recent discussions |2] 

as "superfluid velocities", and that are constructed by dividing the corresponding 3- 

momentum by the relevant mass parameter. The example that arises in the present 

context is that of the so called "superfluid neutron velocity" Vf which is obtained by 

setting 

p\ = mVf. (2.66) 

The corresponding "normal velocity" V 1 is not a pseudo velocity but a true velocity 
in the sense of being the space projected part of a physically well defined 4-velocity, 
namely that of the crust reference frame as given by , so that we can make the 
identification V % = V % c . This leads to a corresponding density decomposition 

P = P s + Pxi ( 2 - 67 ) 

(which unlike (|2.57|) is not affected by the choice of the gauge parameter a c as will 
be shown later) in terms of which the dynamical energy will be expressible simply as 

tfdyn = \{Ps + Pn O , (2-68) 

with the so called "superfluid density" p s and the so called "normal density" p N 
given in the notation of (|2.65J) by 

p s = p { m/m i k , p N = p c + p f m f c /m*. (2.69) 



3 Underlying physical theory of neutron conduction 

3.1 Microscopic analysis in the crust frame 

As it is not practical to carry out direct experimental measurements for bulk matter at 
neutron star densities, the evaluation of the relevant equation of state functions requires 
a theoretical analysis based on an underlying physical model. In particular, to provide 
the quantitative estimates ^Hj needed for obtaining the entrainment density, or the 
corresponding neutronic effective mass in the ionic crust layers above the neutron drip 
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density threshold, it is necessary to use a microscopic model of the kind whose essential 
elements were presented in our preceding work ^H] and that we have recently developed 
in a follow up JH] allowing for the effect of BCS pairing. 

The approach we use provides a typical continuum fluid description as formulated in 
terms of mesoscopically homogeneous configurations that are obtainable from a corre- 
sponding microscopic theory by the minimisation subject to relevant constraints of the 
energy density U that is defined in terms of a quantum system in a unit volume sample 
box (subject to periodic boundary conditions) as the averaged expectation value, which 

we indicate by angle brackets, of the relevant total Hamiltonian operator H say, i.e. 

U = (Hj . The analysis will be carried out using the crust rest frame, in which the 
heavy ionic nuclei can be treated as classical particles at fixed positions, to which all the 
protons and a subset of "confined" neutrons are considered to be bound. For simplic- 
ity a quantum description is applied only to "free" (unconfined) neutrons (whose role 
is analoguous to that of "conduction" electrons in ordinary solid state theory) whose 
Bragg scattering by the potential wells associated with the nuclear clusters gives rise 
to strong entrainment effects. 

In a strictly static configuration the only relevant constraint is the preservation of the 

number density n of such unconfined neutrons, which is definable as the corresponding 

average n = (n) of a conserved particle number operator fl say. The minimisation 

of U subject to such a constraint is equivalent to the absolute minimisation of a 

combination of the form 

U' = U-fin, (3.1) 

in which fl is a Lagrange multiplier that will be interpretable as the relevant chemical 
potential. For any given value of n the corresponding constrained minimum state can 
be expected on symmetry grounds to be static, in the sense of containing no relatively 
moving currents. To obtain the relatively conducting configurations in which we are 
interested here, we need to impose the further constraint of preservation of the relevant 
current components ri 1 (with respect to the rest frame of the sample box, i. e. in the 
crust frame) as defined by the corresponding averages of the relevant quantum operators 
n l = (fl 1 ) say. The problem will thus be that of minimising a combination of the 
form 

U' {p} = U'- Pl rf (3.2) 

in which the quantities Pi are Lagrange multipliers that will be seen to be interpretable 
as representing effective momentum per particle, so that an infinitesimal variation in 
the neighbourhood of such a conditionally minimised energy density will be given by 

SU = nSn+piSn*. (3.3) 
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This means that the multipliers jl and pi can be construed as partial derivatives 
with respect to number density and current respectively. 

Subject to the realistic assumption that the current is small, the energy density will 

be close to the static internal energy value U i ns say (minimising U ) in the state 

of zero current characterised by a given value of the number density n and the ionic 

charge number Z . It will therefore be expressible by an expansion that will be given 

to second order in the current density n l by a decomposition of the form 

U = U ins + U dyn , (3.4) 

in which - for a given value of the charge number Z and the confined particle number 
density n c which, in the present case are considered as given classical variables in the 
microscopic model - the static internal energy density U[ ns will depend just on the 
relevant particle number density n , so that it will be expressible as 

Ui ns = U ins {Z, n c , n} , (3.5) 

while the dynamical contribution is given in terms of a positive definite matrix with 
tensor components /C that also depend just on the same three scalar variables by an 
expression having the homogeneous quadratic form 

^dyn = |/C7/nV . (3.6) 

By differentiating this expression with respect to the current components, it can be seen 
that in terms of the momentum components appearing as partial derivative coefficients 
in the variation formula (jH.Hjl , the current will be given by the formula 

n l = K ij pj . (3.7) 

In practice n and n c can be expected to remain close to values that are determined 
as functions of Z and the total baryon number density n\, by a condition of chemical 
equilibrium whereby U i ns is minimised, with value £7 eq say, for the given values of 
Z and Tib , where n will have a value, n e q say, that is determined as a function 
of Z and 77-b and hence of the corresponding value of n c . Therefore rather than 
considering the coefficients JC 1 ^ to be functions of the three independent variables Z 
n c and n, it should be an adequate approximation to consider them to be functions 
just of the pair of variables Z and n c that determine the classical background field in 
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the model. This allows us to deduce from (J3.3j) that the partial derivative with respect 
to n of the static internal energy density function (j3.5j) will be given exactly by 

-^ = m- (3.8) 

To accuracy of linear order in deviations from chemical equilibrium it will therefore 
suffice to take the function ()3.5|) to have the simple form 

Ums = U eq + fj,(n - n eq ) (3.9) 

in which U e q , 71 e q and jl itself are all functions just of the pair of variables n c 
and Z that characterise the mean density and degree of clustering of the underlying 
distribution of protons. 

The mobility tensor }C lJ defined by the formula (j3.7j) might be anisotropic in a 

solid [Hj of perfectly regular crystalline type, but one would expect it to be given by 

an expression of the isotropic form 

JC tJ = JCY J , (3.10) 



in which the scalar coefficient will evidently be given by 

K = \^V\ (3.11) 



not only for a medium that is a liquid (as will be the case in a neutron star crust 
when the star is very young) and for the case of a solid having a glasslike or disordered 
crystalline structure on a macroscopic scale (as is likely to be the case in a realistic de- 
scription of a neutron star crust) but even for a perfectly regular crystal lattice provided 
it is of a cubic type (such as will be energetically favoured throughout the neutron star 
crust except near the base where rod or slab like structures may be preferred). 

In terms of the mobility scalar introduced in this way, the relation between current 

and momentum (in the crust frame) will be given by 

n* = JCf j pj, (3.12) 

and the corresponding final result for the dynamical energy contribution in (jH.fi j) will 
be given by 

U d yn = ^^PiPj = 2^"^' n%TV> ' (3.13) 
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3.2 Matching with the macroscopic description 

In order to match the microscopic quantum mechanical description in the previous 
subsection to the macroscopic fluid description summarised before, an obvious first 
step is to take the convention that the fraction of the neutrons that is considered to 
be bound to the ionic nuclei is such that we can identify the unbound neutron nunber 
density n with the free neutron number density n introduced above. 
Having thus set 

n f = n (3.14) 

we go on to observe that in the crust frame characterised by V* = in which the 
preceding microscopic analysis was carried out the free neutron current will just be the 
same as the (gauge independent) conduction current ()2.53|) so we shall be able to make 
the further identification n l = Tl\ . 

It can be seen that with respect to this particular frame the kinetic energy will 
simply be given by 

Ukm = lmn { jijVfv£ = \ — ^ lJ n l n J . (3.15) 

77-f 

In order to evaluate the (frame independent) entrainment energy that is definable, in 
accordance with ()2.58|) as 

U en t = U d yn ~ U kin , (3. 16) 

it now suffices to substitute the formula ()3.6|) for the dynamical part, which provides 
the required result in the form 

^e n t = i(^/C7/-^ 7 ^nV, (3.17) 

which is noteworthy for remaining valid for application to the elastic solid models ^3] 
that will ultimately be needed, in which the mobility tensor /C ZJ need not always be 
isotropic. 

In the fluid case characterised by (|3.10J) with which we are concerned here, the 
formula (j3.17|) reduces to 

By comparison with ()2.60|) the entrainment density can be read out from this as 

Pic = nf 2 /JC - pi , (3.19) 
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which by (J2.61|) a implies 



(3.20) 



This quantity depends on the choice of chemical basis, but we can use ()2.56|) to obtain 
the corresponding coefficient 

/C ff = /C~ 1 , (3.21) 

which has the important property of being gauge invariant in the sense that it is 
independent of the choice of the gauge parameter a c that has been used in ()2.23|) for 
specifying the number density Tl{ of neutrons that are counted as "free" . The gauge 
independence of the tensor fC i • and thus of the scalar JC is an evident consequence 
of the gauge independence of the relative current vector ()2.54j) in terms of which it was 
defined by ()3.6|) . 

We can go on to rewrite the formulae ()2.(i9j) for the so called "superfluid" and 
"normal" density contributions in the manifestly gauge independent form 

p s = p — p N = m 2 /C . (3.22) 

This shows that the "superfluid" density contribution p s will tend to zero when 
the mobility coefficient /C tends to zero, which is what occurs at the "neutron drip" 
transition (whose role is thus analogous to that of the "lambda point" transition in 
ordinary liquid helium). 

The actual evaluation of the required mobility coefficient KL as a function of the 
relevant densities has been initiated in the preceding work [H)] which describes the way 
in which a simple nuclear physical treatment can be applied most easily to the rod 
or plate type configurations that are likely to be relevant at the base of the crust [23 
El], while more recent work using more elaborate numerical analysis has extended the 
range of this treatment to include three dimensional cubic configurations as well [TH] . 
The quantitative results obtained so far are of an approximate provisional nature, and 
much further work will be needed to improve their precision and reliability. A first 
step towards such refinement has been taken in our recent examination ^H] of the - 
apparently only moderate - adjustment of the mobility tensor that is needed to take 
account of the pairing effect responsible for the BCS type superfluidity of the neutrons 
that is predicted to occur in the low to moderate temperature range that is relevant in 
pulsars. 

When the appropriate value of /C has been obtained, the complete determination 
of the model still requires the specification of the (frame independent) internal static 
contribution ()2.4()j) . This can evidently be obtained just by adopting the prescription 

_____ r 

(13.9)1 , which enables us to make the simple identification p o — > — p in the zero 
current limit. 
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4 Restoration of frame covariance 



4.1 Chemically invariant action formula 

In the crust rest frame used in the previous section the "free" neutron momentum 
covector was just the quantity P{ that is given by (|3.12j) , so with respect to a generic 
frame with non vanishing 3 velocity V l c this quantity will be given by the formula 

/J. = mjijvj' + pi , (4.1) 

in which it is to be observed that both terms are independent of the parameter a c that 
specifies the chemical gauge, so that as remarked in Subsection 12.21 their is no ambigu- 
ity in the application of the superfluidity condition to the effect that the momentum 
covector (j4.1j) should be irrotational. 

This contrasts with the chemical gauge dependent status of the corresponding 

"crust" momentum covector which will be expressible in terms of the (chemically gauge 

invariant) relative conduction current Tl l (J2.53)) by 

= m lij v i + — Hj n Pi ■ ( 4 - 2 ) 

in which the relative conduction current n l and the corresponding momentum con- 
tribution Pi are gauge invariant, but the densities Tif and n c are not. 

Similar remarks apply to the complete dynamical action density (|2.41|) which can 

be given by the formula 

U dyn = \p Vc + /C (m V l c pi + \p 2 ^j , (4.3) 

in which each separate term is chemically invariant, in contrast with the status of the 
separate kinetic and entrainment contributions in the decompositions ()2.58|) , and in 
the corresponding decomposition 

Amat = Aint + A kin (4.4) 

of the material action contribution in which the kinetic part is simply given by Aki n = 
U kin and the internal contribution will take the form 

A int = f/ent-^ins- (4.5) 

This quantity needs to be distinguished from the corresponding internal energy density 
defined by 

^int = CAnat ~ ^kin , (4.6) 
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which will be given by 

^int = ^ent + ^ins. (4.7) 

In particular it can be seen from (|4.5jl and ()4.7|) that the internal contribution Aj n t 
to the Lagrangian density is not simply given by the negative of the internal energy 
density Ui n t (as stated by Prix et al. El El El) but takes the form 

Aint = -£Ant + 2?7ent • (4.8) 

The tilde is used here as a reminder that the contributions thus designated have the 
disadvantage of being chemically gauge dependent. However the terms in (14. 5 J) have the 
compensating advantage of Galilean frame invariance. It is therefore convenient and 
customary to use these entities as a starting point for the specification of particular 
models, as they contain just the minimum information, in the form of the particular 
equations of state, that are needed for this purpose, namely expressions as a function 
of the relevant scalar densities Tlf , n c and if necessary also ll\ , for the static internal 
energy U j ns and for the mobility coefficient JC or equivalently, via (j3.19j) , of the 
entrainment density p^ c for the entrainment part U en t . In terms of the relative 
velocity defined by ()2.52|) the generic variation of the frame invariant internal action 
density will evidently be expressible in the form 

£Aint = p ic V{ c i Svf c - x f £n f - x c ^ c - ffini . (4.9) 

The frame independent chemical potentials X X ■> X defined in this way determine 
corresponding potentials 

M f = X f - \m v { 2 , ji c = x° ~\m vl , Ji l = x 1 , (4. 10) 

(of which the first two are Galilean frame dependent) in terms of which the variation 
of the complete (chemical gauge invariant) material action density (j4.4j) will take the 
form 

^A mat = riffJ^v^ + ric/i^Vc — J^5n{ — jl c 5n c — p^Sni . (4-11) 

To match this with the preceding expression ()2.24j) for the material action variation, 
it can be seen that we need to make the identifications 

M f = -UfV„ , A c = , (4.12) 
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(whereas to relate this to the terminology of Prix et al. [S] one would need to make 
the translations X — > X° — > an d -A-int — > —£ )• 

To match this with an alternative treatment of a kind suitable for generalisation 
[Tl] to allow for solid elasticity, in which the confined constituent will have a privileged 
role, the relative velocity needs to be replaced in the variation (|4.9J1 by the relative 
current n l given by ()2.53|) which has the advantage of being chemically invariant. In 
the notation of (|2.65j) one thus obtains 

f 

M int = — - lij n l 5n J - xSn { - x c $n c - x^i , (4.13) 

with 

X f = % f + m f c ^ 5 x c = X c , X l = X l - (4-14) 
4.2 Effective mass relations 

Subject to the understanding that number density n in the preceeding work 

is to be interpreted as the "free" number density 7lf in the present terminology, the 

corresponding effective mass will be identifiable with the effective mass per free 

particle as denoted here by , which means that it will be given according to (j3.21j) 

by 

ml = n { /K. (4.15) 

A particularly noteworthy and previously unexpected conclusion that has emerged 
from our preliminary investigations EI] is this effective mass can be expected to 
become very large (by a factor of several hundred percent) compared with the ordinary 
neutron mass m in the middle layers of the outer crust, even when it is only the 
neutrons outside the nuclei that are counted as free. It would of course be even larger 
with respect to the simple "comprehensive" gauge given by a c = 1 for which all the 
neutrons are counted as "free". 

Before concluding, it needs to be pointed out that whereas the definitions used here 

for the effective masses are equivalent to setting 

ml = JC s n { , ml = /C cc n c , (4. 16) 

there is an alternative definition that is commonly used in the published literature on 
neutron star matter [3]. This alternative definition can be formulated in the present 
terminology by setting 

m\ = r?f / JCs , m\ = nj JC CC . (4. 17) 
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The effective masses of this second kind will be given in terms of those of the first kind 
by 

raj — m = — - (m^ — ra) , m| — m = — j (m° — ra) . (4. 18) 

These effective masses of the second kind are (like the first kind) either both larger or 
both smaller than the ordinary mass. 

To complete the specification of the two fluid description, the coefficient K} c giving 

rise to entrainment can be seen to be given by either of the equivalent formulae 

K c n c = — ra f c = ra — , K c ri{ = —m c { = m — m£ , (4.19) 

Thus K^ c will vanish, as expected, whenever the effective masses are equal to the 
ordinary (protonic) baryon mass ra . (As previously remarked, both effective masses 
are equal to ra whenever one of them is equal to ra ). The alternative description in 
terms of the effective masses of the second kind leads to the expressions 

m ml m 

The conditions ()2.49|) can thereby be restated in terms of effective masses in the 
equivalent forms 

ml n c ml m , A 

— >— , — >— , 4.21 

m rz-b m nb 

while the effective masses of the second kind must obey the equivalent inequalities 

m^nj, S < ^. (4 22) 

m n c m ri{ 



4.3 Comprehensive treatment of inner crust and (outer) core 

The approach developed above has been particularly designed for the treatment of the 
inner crust layers above the "neutron drip" density threshold but below the transition 
density at which the ionic clusters merge. It is desirable to relate this approach to that 
of related work0 El \7\ designed more specifically for treatment of the homogeneous 
(unclustered) nuclear matter of the core region, or to be more precise the outer core 
region (since the most massive neutron stars may also contain an inner core consisting 
of hyperons or quark matter that would need a relativistic treatment of kind that is 
beyond the scope of the present discussion). 
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Liquid core 



Inner crust 



Figure 1: Rough sketch of expected radial dependence of dimensionless "comprehensive" 
entrainment coefficients £ n (dashed line) and £ p (solid line) in a neutron star. 

Instead of using the more general kind of chemical basis that is useful for discussion 
of the "neutron drip" transition, for the purpose of matching to the formalism that 
is usually employed for work on the outer core it is convenient to work exclusively in 
the "comprehensive" gauge characterised by setting CL C = 1 in (|2.23j) , which simply 
means identifying the "free" neutron current 71^ with the entire neutron current 
and identifying the "confined" baryon current 71^ with the proton current Tip v . 
This allows us to obtain a treatment combining the representation of the crust and the 
core in a single "comprehensive" model in which the relevant "comprehensive" effective 
mass, 

ml = njK, (4.23) 

(of the first kind) for the neutrons has the advantage of having a clear physical def- 
inition, though - as the price for this - it has the awkward property of diverging to 
infinity as the density decreases to that of the "neutron drip" threshold. 

The link with work on the core regime can be made more explicit by using notation 

of the kind employed by Prix et al. [HI El Hj in which the neutron and proton 3- 

momenta are given in terms of the relative 3-velocity Z\ z = V 7 ^ — V l n by expressions 

of the form 

iA = nM v i + £ n Aj ) > A = ^M v i - ^ AJ ) > ( 4 - 24 ) 

involving dimensionless parameters £ n and £ p that can be seen from ()2.44j) and 
(jPU) to be definable as 

e n = (m — m+)/m , e v = (m — m£) / 'm . (4.25) 
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Liquid core 



Inner crust 



Figure 2: Rough sketch of expected radial dependence of dimensionless "operational" 
entrainment coefficients £f (dashed line) and £ c (solid line) in a neutron star. 

According to an identity of the same kind as that given by (|4.19jl , these coefficients 
S n and £ p will be related to each other, and to a density variable a introduced by 
Prix et al. by 

£nn n = £ P n p = 2a/m , a =-\p np . (4.26) 

This means in particular that for any given layer inside the star, the coefficients £ p 
and £ n both have the same sign, and it also shows that we shall have |e n | < |e p | 
since the density of neutrons will always exceed the number of protons, increasing from 
a comparable value in the outer layers, where we shall have £ n ~ £ p , to a very much 
higher value, n n ^> n p in the neutron rich core where we shall have \s n \ |£ p |- In 
terms of the proton fraction £ p = Tip/ 71^ < 1/2, it can be seen that an inequality 
of the type (|4.21j) is equivalent to either of the restrictions 

e n < x p , e p < 1 - x p . (4.27) 

which evidently entails that the entrainment coefficients £ n and £ p must both be less 
than unity. 

As the effective mass (and hence also ) has been found [2211201 to be smaller 
than the ordinary nucleon mass 772 in the outer core, it follows that the entrainment 
coefficients £ n and £ p will both be positive there. This contrasts with situation 
in the crust, where we expect the "comprehensive" effective mass ?7i£ > m to be 
larger compared with m so that the coefficients £ n and £ p will both be negative, 
since we have found [TI)l ITU] that the smaller effective mass given in a "realistic" gauge 
just for "conduction" neutrons with number density rif (defined as those contributing 
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Figure 3: Rough sketch of expected radial dependence of "comprehensive" effective 
mass ratios m^/m (dashed line) and m*/m (solid line) in a neutron star. 

significantly to the current) by = Tn^n{/n n will already be larger (and in some 
layers much larger) than the ordinary mass TTl in the crust. It is to be observed that 
the negative values of £ n and £ p will actually diverge at the "neutron drip" transition 
where the "conduction" number density 7if tends to zero. These considerations are 
summarised on figures d El 01 and |U 

The total material energy density (as locally defined without allowance for gravity) 

will be expressible in the "comprehensive" treatment by 

U ma t = C^kin + U' m t , U m t = U[ ns + U en t > (4.28) 

with the kinetic contribution given by the usual formula 

^kin = ^m(n n vl + n v v v 2 ) , (4.29) 

while from equations ()2.60|) , ()2.65|) and (|4.25J) it can be seen that the two fluid equa- 
tion of state will be expressible in terms of a pair of 3- variable functions U m s{^I? n n, n p} 

and Oi{n\, n n , n p } in terms of the quantity that has been loosely referred by Prix et 
al. [3 El EE Hj to as the "internal energy density" - but that is actually the negative of 
the internal action density - will be given by 

-Aint = ^ins + «^ 2 . (4.30) 

whereas the true internal energy density will be given by 

U int = U- ms -aA 2 . (4.31) 



28 



o 1 

Liquid core 



1 — ► 

Inner crust 



Figure 4: Rough sketch of expected radial dependence of "operational" effective mass 
ratios, m^/m (dashed line) and m^/m (solid line) in a neutron star. 

In the high density core region the static energy contribution U ins will cease to have 
any dependence on the ionic number density n\ , not because of any tendency of Tl\ to 
vanish on the lower boundary of the crust (where it merely ceases to be well defined) 
but rather because the corresponding clustering energy coefficient jl will vanish there, 
so that the value of TL\ no longer matters. 

As depicted in figure ^ the entrainment coefficients between the neutron and pro- 
ton fluids in the core of a neutron star are expected to be very small compared to 
the previously neglected entrainment coefficients in the inner crust layers. This strong 
entrainment in the inner crust may have significant effects on the dynamical evolution 
of the star. For instance, as recently shown by Andersson et al. sufficiently large 
negative values of the entrainment coefficients £ n and £ p could trigger a Kelvin- 
Helmholtz instability in the fluid mixture which might be at the origin of observed 
pulsar glitches. The development of gravitational wave detectors could provide valu- 
able constraints upon the internal structure of a neutron star by asteroseismology |2Z| • 
However a thorough understanding of the oscillation modes of a neutron star is required 
in order to extract useful information from the gravitational wave data. Hydrodynam- 
ical simulations of neutron star cores indicate that these oscillation modes are very 
sensitive to entrainment [22] . The inner crust, where the two fluids are expected to 
be the most strongly coupled, should therefore be given careful consideration in neu- 
tron star dynamical studies. As shown in the present work, a simple neutron star two 
fluid model allowing for the presence of the crust (but neglecting the stress anisotropy), 
can be easily implemented by suitable adjustments of the corresponding entrainment 
coefficients. 
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